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ABSTRACT 

Recent proper motion and parallax measurements for the pulsar PSR B1508-I-55 in- 
dicate a transverse velocity of 1 100 km s^^, which exceeds earlier measurements for 
any neutron star. The spin-down characteristics of PSR B1508-I-55 are typical for a 
non-recycled pulsar, which implies that the velocity of the pulsar cannot have orig- 
inated from the second supernova disruption of a massive binary system. The high 
velocity of PSR B1508+55 can be accounted for by assuming that it received a kick 
at birth or that the neutron star was accelerated after its formation in the supernova 
explosion. We propose an explanation for the origin of hyperfast neutron stars based 
on the hypothesis that they could be the remnants of a symmetric supernova explo- 
sion of a high- velocity massive star which attained its peculiar velocity (similar to that 
of the pulsar) in the course of a strong dynamical three- or four-body encounter in 
the core of dense young star cluster. To check this hypothesis we investigated three 
dynamical processes involving close encounters between: (i) two hard massive bina- 
ries, (ii) a hard binary and an intermediate-mass black hole, and (iii) a single star 
and a hard binary intermediate-mass black hole. We find that main-sequence 0-type 
stars cannot be ejected from young massive star clusters with peculiar velocities high 
enough to explain the origin of hyperfast neutron stars, but lower mass main-sequence 
stars or the stripped helium cores of massive stars could be accelerated to hyperveloc- 
ities. Our explanation for the origin of hyperfast pulsars requires a very dense stellar 
environment of the order of 10® — 10'' starspc"'^). Although such high densities may 
exist during the core collapse of young massive star clusters, we caution that they 
have never been observed. 

Key words: stellar dynamics - methods: N-body simulations - stars: individual: RX 
J0822-4300 - stars: neutron - pulsars: general - pulsars: individual: B1508-f 55. 



1 INTRODUCTION 

It has been known for a long time that the typical pecu- 
liar space velocities of pulsars are an order of magnitude 
larger that those of their progenitors - the massive stars 
(e.g. Gott, Gunn & Ostriker 1970). Subsequent proper mo- 
tion measurements for pulsars suggested that some of them 
could be very fast objects, with peculiar speeds of up to 
~ lOOOkms"^ (e.g. Chatterjee & Cordes 2004; Hobbs et al. 
2005). Two mechanisms have been proposed for the origin of 
high- velocity pulsars. The first one involves the disruption 
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of tight (semidetached) massive binary systems following 
the second (symmetric) supernova explosion (e.g. Blaauw 
1961; Iben & Tutukov 1996). This mechanism, however, can- 
not produce velocities in excess of ~ 700 kms"^ (Porte gies 
Zwart & van den Heuvel 1999). Moreover, the high- velocity 
pulsars formed in this way are the remnants of the first su- 
pernova explosion and therefore should be recycled, i.e. their 
spin characteristics should be significantly affected by the 
stellar wind of the companion star - the progenitor of the 
second supernova (cf. Chatterjee et al. 2005). The second 
mechanism relies on a natal kick (caused by the asymmetry 
of the supernova explosion; e.g. Shklovskii 1970; Dewey & 
Cordes 1987) or on a post-natal acceleration (caused by the 
asymmetric electromagnetic or neutrino emission from the 
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newborn pulsar; e.g. Tademaru & Harrison 1975; Chugai 
1984). In principle, this mechanism can produce solitary 
high-velocity pulsars with ordinary spin-down characteris- 
tics, i.e. non-recycled pulsars (e.g. Scheck et al. 2006). Thus, 
detection of non-recycled pulsars moving with a velocity 
^ lOOOkms"^ would suggest that the mechanism involv- 
ing the binary disruption by the supernova explosion does 
not apply to all cases. In the following we refer to pulsars 
(neutron stars) moving with velocities of > 1000 km s^^ as 
'hyperfast pulsars' (a term introduced by Chatterjee et al. 
2005). 

Until recently the best candidates for hyperfast pulsars 
were PSR B2224-I-65 (associated with the well-known Guitar 
Nebula; Chatterjee & Cordes 2004) and PSR B2011-h38 (see 
Hobbs et al. 2005 and references therein) , whose peculiar ve- 
locities (both ~ 1500kms~^) were inferred on the basis of 
proper motion measurements and dispersion measure dis- 
tance estimates. The uncertainties associated with the dis- 
tance estimates, however, leave a possibility that these veloc- 
ities are overestimated. Recent proper motion and parallax 
measurements for the pulsar PSR B 1508-1-55 gave the first 
example of a high velocity (1 OSS^gg"^ kms^^) directly mea- 
sured for a (non-recycled) pulsar (Chatterjee et al. 2005). 
This result proved the existence of a population of hyper- 
fast neutron stars, and implies that the peculiar velocity of 
PSR B 1508-1-55 cannot be solely due to the disruption of 
a tight massive binary system. A possible way to account 
for the high velocity is to assume that at least part of this 
velocity is due to a natal kick or a post-natal acceleration 
(Chatterjee et al. 2005). 

In this paper, we propose an alternative explanation for 
the origin of hyperfast pulsars (cf. Gvaramadze 2007). We 
suggest that PSR B1508-I-55 (as well as other hyperfast neu- 
tron star^ could be the remnant of a symmetric supernova 
explosion of a high- velocity massive star (or its helium core) 
which attained its peculiar velocity (similar to that of the 
pulsar) in the course of a strong three- or four-body dynam- 
ical encounter in the core of the parent young massive star 
cluster (YMSC). 

Our suggestion is based on the fact that massive (^ 
10* M0) star clusters are still forming in the disk of the 
Milky Way (see Sect. [3} and on the hypothesis that the 
cores of some YMSCs harbour intermediate-mass black holes 
(IMBHs), i.e. black holes (BHs) with masses ranging from 
~ 100 to ~ 10* Mq (see Sect.[21). Our suggestion also re- 
quires the existence of very dense stellar environment (~ 
10^ — lO'^ pc~^) to ensure that three- or four-body dynami- 
cal encounters are frequent (see Sect.U). Although such high 
densities may exist during the core collapse of YMSCs, we 
caution that they have never been observed (we discuss this 
problem in Sect.|3J. In Sect. [51 we briefly discuss the mecha- 
nisms for the origin of hypervelocity stars, the recently dis- 
covered class of stars moving with a speed of ~ 1 OOOkms"^, 
and suggest that similar mechanisms (although acting in a 

^ Another example of a hyperfast neutron star was recently re- 
ported by Hui & Becker (2006) and Winkler & Petre (2007). Their 
proper motion measurements for the central compact object RX 
J0822-4300 in the supernova remnant (SNR) PuppisA suggest 
that the peculiar velocity of this neutron star could be as large 
as ~ 1000 — 1500kms~^, provided that the distance to these 
objects is ~ 2 kpc. 



different environment) could be responsible for the origin 
of some of the hyperfast neutron stars. In Sect.|4l we esti- 
mate the maximum possible velocities of high- velocity esca- 
pers produced by various dynamical processes in the cores of 
YMSCs. In Sect.O we compare these estimates with results 
from numerical simulations while in Sect.|S] we discuss the 
ability of these processes to explain the origin of hyperfast 
pulsars. 



2 HYPERVELOCITY STARS AND 

MECHANISMS FOR THEIR PRODUCTION 

The discovery of high velocity pulsar PSR B1508-I-55 (Chat- 
terjee et al. 2005) coincided with the discovery of the hy- 
pervelocity starqj (Brown et al. 2005; Edelmann et al 2005; 
Hirsch et al. 2005). The existence of the latter was predicted 
by Hills (1988), who showed that a close encounter between 
a tight binary system and the supermassive BH in the Galac- 
tic centre (e.g. Chez et al. 2003; Schodel et al. 2003) could be 
responsible for the ejection of one of the binary components 
with a velocity of up to several lOOOkms"'^ (see also Hills 
1991). Yu & Tremaine (2003) proposed two additional pos- 
sible mechanisms for the production of hypervelocity stars. 
The first one involves close encounters between two single 
stars in the vicinity of a supermassive BH (the probability of 
this process, however, is very low and we will not discuss it 
further) while the second one is based on the interaction be- 
tween a single star and a putative binary BH in the Galactic 
centre (e.g. Hansen & Milosavljevic 2003). 

The discovery of hypervelocity stars (Brown et al. 2005; 
Edelmann et al. 2005; Hirsch et al. 2005; Brown et al. 2006) 
provides support for models of ejection mechanisms involv- 
ing dynamical processes in the vicinity of the supermassive 
BH (Gualandris, Portegies Zwart & Sipior 2005; Baumgardt, 
Gualandris & Portegies Zwart 2006; Bromley et al. 2006), 
and now it is widely believed that these high-velocity objects 
originate in the Galactic Centre (see, however, Edelmann et 
al. 2005; see also Sect.[6}. It is, therefore, possible that some 
hyperfast pulsars (or their progenitor stars) were also ejected 
from the Galactic Centre. For example, one of the hyperve- 
locity stars is a early B-type star of mass ~ 8 Mq (Edelmann 
et al. 2005) , so that this star will end its evolution in a Type 
II supernova leading to the production of a hyperfast neu- 
tron star. The proper motion measured for PSR B1508-(-55, 
however, indicates that this pulsar was born in the Galactic 
disk near the Cyg OB associations (Chatterjee et al. 2005), 
i.e. its origin cannot be associated with the Galactic Centrqj. 
The proper motions of PSR B2224-I-65 and PSR B201H-38 
also suggest that these pulsars were born far from the Galac- 
tic Centre. Therefore, one should look for other places in 
our Galaxy where one can find a sufhciently massive BH 
and where the number density of the local stellar popula- 
tion is high enough to ensure that close encounters between 
stars and the BH are frequent. Note that the characteristic 

^ In the following we use the term 'hypervelocity' [introduced by 
Hills (1988)] to designate the ordinary stars moving with a pe- 
culiar speed of > 1 000kms~^ while for neutron stars we reserve 
the term 'hyperfast' (see Sect. 1). 

For a nice picture illustrating the trajectory of the pulsar on the 
sky see |http : //www. jb. man. ac.uk/news/f astestpulsar.html | 



age of the above mentioned pulsars (ranging from ~ 0.4 to 
~ 2.3 Myr) implies that at least these three pulsars were 
not ejected from the dense cores of globular clusters (whose 
typical age is > Gyr) and point to the more plausible cites 
of their origin - the YMSCs. 

We note that with the YMSCs could be associated 
another important channel for production of high-velocity 
stars, namely through close dynamical encounters between 
stars in their cores. This process constitutes the base of the 
dynamical-ejection scenario proposed by Poveda, Ruiz & 
Allen (1967) to explain the origin of runaway OB stars. The 
most effective path for production of high- velocity stars by 
stellar encounters is through interaction between two hard 
binaries (e.g. Mikkola 1983; Leonard & Duncan 1988, 1990). 
This process and processes involving dynamical encounters 
with IMBHs are discussed in detail in Sect.U] 



3 YOUNG MASSIVE STAR CLUSTERS 

Recent discovery of young (^C lO'^ yr) and massive 
10* Mq) star clusters (Borissova et al. 2006; Figer et al. 
2006; Davies et al. 2007) accompanied by upward revision of 
masses of the already known star clusters (Knodlseder 2000; 
Clark et al. 2005; Homeier & Alves 2005; Santos, Bonatto & 
Bica 2005; Ascenso et al. 2007; Wolff et al. 2007; Harayama, 
Eisenhauer & Martins 2008) increased the number of YM- 
SCs in the Galactic disk to > 10. All but one of these YMSCs 
are located on the near side of the Galaxy. Taken together 
these facts suggest that YMSCs are more numerous than it 
was known hitherto and that many of them are still hidden 
from observers by the obscuring material in the Galactic 
plane (Knodlseder et al. 2002; Hanson 2003). 

Let us estimate the number of YMSCs formed in the 
Galactic disk during the last 3 x 10^ yr (the lifetime of a 
star with mass of 8 Mq - the minimum mass of single stars 
producing neutron stars) . We choose this time-span since we 
assume that the currently observed hyperfast neutron stars 
could be the descendants of either a 8 M0 hypervelocity 
star, ejected at the very beginning of its life, or a stripped he- 
lium core of the more massive star, attained its peculiar ve- 
locity only recently (see Sect.Q. Assuming that the star for- 
mation rate (SFR) of ~ 7-10x10"* M© yr"^kpc"^ derived 
by Lada & Lada (2003) for embedded star clusters in the so- 
lar neighbourhood, is representative for the Galactic disk as 
a whole (cf. Larsen 2006), one has that ~ 6.6-9.4X lO'' Mq 
were formed within a circle of 10 kpc. For a power-law clus- 
ter initial mass function with a slope —2 (e.g. Lada & Lada 
2003) and 50 Mq < M^ < 10** Mq, one has ~ 70-100 YM- 
SCs with a mass Md ^ 10* Mq . This estimate is in good 
agreement with that derived by Knodlseder et al. (2002) us- 
ing different arguments (see also Madhusudhan et al. 2007) , 
and exceeds by a factor of 2-3 the figures derived by Larsen 
(2006). The discrepancy with Larsen (2006) originates due 
to the somewhat different input parameters adopted in his 
paper. For example, he uses a factor of 1.3-1.9 smaller SFR, 
derived by Lamers et al. (2006) for bound star clusters in the 
solar neighbourhood, and assumes a larger upper cut-off of 
cluster masses, Md < 10^ Mq. We prefere to use the SFR 
given in Lada & Lada (2003) to account for the fact that 
many embedded star clusters independently of their mass 
became dissolved within the first ^ 10^ yr of their evolution 
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(a possible example of such a stellar system is the Cyg OB2 
association) . 

To the YMSCs resided in the disk one should add the 
ones formed in the central part of the Galaxy. Theoretical 
models of evolution and observability of star clusters within 
200 pc of the Galactic Centre suggest that this region could 
easily harbour ~ 50 YMSCs with properties similar to those 
of the well-known Arches and Quintuplet systems (Portegies 
Zwart et al. 2001a, 2002). Recent discovery with 2MASS 
of a large number of cluster candidates in direction to the 
Galactic Centre (Dutra et al. 2003, and references therein) 
supports this suggestion, although follow-up work is required 
to confirm the nature of these objects. 

The dynamical evolution of YMSCs is dominated by 
massive stars, which sink to the cluster centre on a time scale 
tec ~ 0.1 — 0.2frh, which typically is smaller than the main- 
sequence lifetime of these stars. Interestingly, the dynamical 
friction time-scale is independent of cluster mass, its size 
and density profile (e.g. Portegies Zwart & McMillan 2002; 
Giirkan, Freitag & Rasio 2004), where t^h is the half-mass 
relaxation time-scale given by (Spitzer 1987) 

* - GV2(m)lnA ' ^ ' 

Here, M^i and rh are the total mass and the characteristic 
(half-mass) radius of the cluster, G is the gravitational con- 
stant, (m) = Mc\/Nci is the mean stellar mass, A'ci is the 
number of stars in the cluster and In A ~ 10 is the Coulomb 
logarithm. For tec = 0.15 and using Eq. one has 




As a result, mass segregation drives YMSCs to core col- 
lapse. The consequential increase in the central density may 
results in runaway stellar collisions during which a massive 
star continues to grow in mass due to repeated bombard- 
ment to a very massive star (VMS) (Portegies Zwart et 
al. 1999; Portegies Zwart & McMillan 2002; Giirkan et al. 
2004), provided that core collapse is completed before the 
most massive stars in the cluster explode as supernovae, i.e. 
tec ;$ 3 X 10® yr. On this time-scale essentially all massive 
(> 20 Mq) stars in YMSCs segregate to the cluster core. 
For Mci = 10* - lO'' Mq and assuming a 0.2 - 120 Mq 
Salpeter initial mass function (in this case (m) ~ 0.7 Mq 
and iVci ^ 1.5 x 10* - 10^), one has from Eq. ^ that the 
runaway process occurs in YMSCs with rh < 0.3 — 0.8 pc. 
The majority of young and massive stars appear to reside in 
dense embedded clusters with a characteristic radius of < 1 
pc, which is independent of mass (Kroupa & Boily 2002). 
It is therefore not unconceivable that an appreciable frac- 
tion of YMSCs evolve through a coUisional stage. The maxi- 
mum observed density in the cores of YMSCs [e.g. NGC 3603 
(Harayama et al. 2007) , Arches (Figer et al. 1999)1 is smaller 
than that required for runaway collisions to occuQ. Whether 

* Note, however, that the Arches cluster is not in a state of core 
collapse, but rather on its way towards this (Portegies Zwart et 
al. 2007). 
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or not the discrepancy between the observations and the 
theory is the result of the small number of known YMSCs 
remains unclear to date. The consequences of the coUisional 
stage are profound, and we assume here that it may exist at 
least in some YMSCs. 

The VMS formed through runaway collisions will ulti- 
mately collapse to a BH. One of the main questions here 
is wheather or not the BH will be massive or not. The 
mass of the IMBH formed in this way could be as large as 
~ 1 000 Mq for Mci ~ 10^ M0 (Portegies Zwart et al. 2004), 
but if the mass loss in the collision product is as high as 
argued in some calculations (Belkus, Van Bever & Vanbev- 
eren, 2007; Yungelson 2007; Yungelson et al 2008) the final 
BH may be < 100 Mq. 

It is also possible that IMBHs of mass ~ 10^ - 10* M© 
are the descendants of old globular clusters disrupted by 
the tidal field of the Galactic disk or that they are formed 
from the core collapse of massive population HI stars. Later 
they could be captured gravitationally by molecular clouds 
or by the already existing young star clusters. In the first 
case, the IMBH initiates star formation in the cloud and, 
therefore, resides at the centre of the newly formed star clus- 
ter (Miller & Hamilton 2002), while in the second case the 
IMBH rapidly sinks to the centre of the cluster due to the 
dynamical friction (Miller & Colbert 2004). 



4 ORIGIN OF HYPERFAST PULSARS 

We now consider the possibility that dynamical processes 
in the cores of YMSCs could be responsible for ejection of 
hypervelocity massive stars or their helium cores, whose sub- 
sequent collapse and (symmetric) supernova explosion result 
in the origin of hyperfast neutron stars (pulsars). In the fol- 
lowing subsections, we estimate the upper limits for the ejec- 
tion speed produced by encounters involving: (i) two hard 
binaries, (ii) a hard binary and an IMBH and (iii) a single 
star and a hard binary IMBH. 



4.1 Binary-binary encounters 

In Sect.|2] we mentioned that the most effective path for 
production of high-velocity stars by stellar encounters is 
through interaction between two hard binaries, either the 
tidally captured or the primordial ones. Note that the tidal 
binaries are very hard by definition since the semimajor axes 
of their circularized orbits are at most a few stellar radii (e.g. 
Lee & Ostriker 1986; McMillan, McDermott & Taam 1987). 
Therefore, it is likely that collisions involving tidal binaries 
play a dominant role in production of the fastest runaway 
stars. It is also important to note that due to the mass segre- 
gation, the cores of YMSCs are over-represented by massive 
stars, so that most of tidally captured binaries should be 
the massive ones. Moreover, one can expect that almost all 
tidal binaries would be massive if the most massive stars are 
formed near the centres of their parent clusters (e.g. Bonnell, 
Bate & Zinnicker 1998; cf. Kroupa 2001). 

Numerical simulations by Leonard & Duncan (1988, 
1990) showed that the typical velocities at infinity (i.e. well 
after ejection) of runaway stars (produced in the course of 
binary-binary collisions) are similar to the orbital velocities 



of the binary components while the velocities of some esca- 
pers can be twice as large. Moreover, the maximum possible 
velocity attained by the lightest member of the binaries in- 
volved in the interaction (e.g. the helium core of a massive 
star or a early-type B star) can be as large as the escape 
velocity from the surface of the most massive star in the 
binaries (Leonard 1991). For the upper main-sequence stars 
with the mass-radius relationship (Habets & Heintze 1981) 
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where tms and rriMS are the stellar radius and mass, the 
maximum possible velocity at infinity of ejected stars is a 
weak function of ttims, V^^ms — 700 kms"^(mMs/ Mq)"'^^ 
and could be as large as ~ 1400 kms^^ (cf. Leonard 1991). 
The ejection velocity could be even larger if the binaries 
involved in the interaction are already evolved through a 
common-envelope phase and consist of two helium cores. In 
this case, 
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where rua and mne are the radius and the mass of a he- 
lium core (Tauris & van den Heuvel 2006), so that V™^;, ~ 
14OOkms"^(mHo/M0)'' "^ - 2300kms"\ 

Scattering experiments by Leonard (1991) showed that 
only a small fraction (less than 1 per cent) of runaway stars 
released in the course of binary-binary collisions can attain 
the maximum possible velocity. So that, to produce at least 
one hypervelocity (massive) star during the first several Myr 
of cluster evolution, the coUisional time-scale for binaries in 
the cluster core should be < 10* yr. The total rate of binary- 
binary encounters in the core of radius Vc is 



(5) 



where A^b — (47r/3)rj^nb, rib is the number density of bina- 
ries, Koi is the relative velocity of approach of the binaries 
at infinity. 
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is the gravitationally focused cross section (Leonard 1989), 
and a is the binary semimajor axis. For the sake of simplicity, 
we consider the equal-mass and equal-energy binaries, and 
assume that one of the binary components is a 40 Mq main- 
sequence star while its companion is a He core of mass of 
5 Mq . From Eqs. ((5} and ((Sjl , one has the mean time between 
binary-binary collisions (cf. Leonard 1989) 



tbb ~ 8 X 10^ yr 
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It follows from Eq. (0 that the coUisional time-scale is short 
enough (i.e. < 10* yr) if more than 10 massive binaries ex- 
ist in the core of radius ~ 0.01 pc. Assuming that almost 
all stars more massive than ~ 20 M© are segregated to the 
cluster core and exchanged in binaries, one has the minimum 
mass of the YMSC of ~ 10* Mr. 
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The recent discovery (Muno et al. 2006) of an anoma- 
lous X-ray pulsar in the (4 ± 1) x 10®-yr old massive (> 
10^ M©) Galactic star cluster Westerlundl (Clark et al. 
2005) implies that the progenitor of this neutron star was a 
star with a zero-age main-sequence mass > 40 Mq (cf. Van- 
beveren, Van Rensbergen & De Loore 1998). Therefore, one 
can expect that the hypervelocity helium cores of massive 
stars ejected from the cores of YMSCs during the first sev- 
eral Myr of cluster evolution will end their lives as hyperfast 
neutron stars (pulsars). 



(looll;) (sor;;) • ^^^^ 

It follows from Eq. (|12|) that during the several Myr of clus- 
ter evolution one can expect at least a dozen of close encoun- 
ters involving binaries with a ^ 30 R©. Here, we assume a 
central density in the cluster which exceeds the observed val- 
ues. Although such high density may exist during core col- 
lapse, we note that the well-studied YMSCs in the Galaxy 
are not in this stage of their evolution. 



4.2 Exchange encounters between binary stars 
and an IMBH 

Let us assume that the core of a YMSC harbours an IMBH 
of mass AfsH = 100 — 1 000 Mq , formed either through a 
runaway sequence of mergers or from the core collapse of a 
massive population HI star (see Sect.|3]). 

A close encounter with the IMBH results in the tidal 
breakup of the binary, after which one of the binary com- 
ponents becomes bound to the IMBH while the second one 
(usually the least massive star) is ejected with a high speed, 
given by (Hills 1988; Yu & Tremaine 2003): 

1/4 

(8) 
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where 
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is the tidal radius and mi and m2 are, respectively, the 
masses of the bound and ejected stars (mi > m2). The 
fastest stars produced in this exchange process come from 
encounters involving the tightest binary systems (e.g. tidal 
binaries), since the tighter the binary the closer it can ap- 
proach to the IMBH before tidal breakup. Combining Eqs. 
® and dH), one has 
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For a binary having a massive main-sequence compo- 
nent of mass mi ^ m2 and a ~ 2.5ri,MS (cf. Staffer, Os- 
triker & Cohn 1987), one has from Eqs. ([10} and Q 
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For i\fBH = 100-1 000 Mq and mi = 40 - 100 Mq , one has 
from Eq. ((TTJ that the low-mass binary component (e.g. the 
helium core of a massive star or a B star) can be ejected 
with a speed Vex. ~ 900 — 1300kms~^. The ejection veloc- 
ity could be somewhat higher if the binary involved in the 
encounter consists of two helium cores (cf. Sect. 4.1). For ex- 
ample, assuming that mi.uc = 10 M© and m2,Hc ~ 5Mq, 
and a ~ 5 — 10 R©, one has Vao — 850 — 1 760kms~^. Note 
that the weak dependence of V^o on Mbh implies that the 
hypervelocity helium cores can be produced by exchange 
encounters with stellar mass (i.e. ~ 20 M©) BHs. 

The mean collision time between binary stars and an 
IMBH is (e.g. Binney & Tremaine 1987) 



fcoii ~ 2 X 10" yr 
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4.3 Encounters between single stars and a binary 
IMBH 

Numerical simulations by Giirkan, Fregeau & Rasio (2006) 
showed that runaway collisions and mergers of massive stars 
in YMSCs with initial binary fraction larger than ~ 10 per 
cent could result in the origin of two VMSs. The subsequent 
supernova explosions of these VMSs produce two IMBH, 
which ultimately exchange into a binary (Giirkan et al. 2006; 
cf. Fregeau et al. 2006). 

The IMBH binary (IMBHB) gradually hardens due to 
the interaction with stars in the cluster's core. When the 
binary separation reduces to a < Oh — G^/Aa^, where 
p = MiM2/{Mi + Ah) and Mi and M2 are the component 
masses of the IMBHB (Mi > M2), most of stars passing in 
the vicinity a) of the IMBHB are expelled from the core 
at high velocity. The average ejection speed attained by the 
escapers is (Yu & Tremaine 2003) 
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Some escapers, however, can reach much higher velocities. 
For an IMBHB with mass ratio of ~ 1 the maximum ejection 
velocity is 
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where Vbin ~ [G{Mi + Ah)/a]^'''^ is the relative velocity of 
the binary components if their orbits are circular (Tutukov 
& Fedorova 2005). It is clear that the smaller a the larger 
(Voo) and VJ^^^". There are, however, two constraints on the 
minimum value of a, which limit the maximum possible ejec- 
tion velocity. 

First, a should be sufficiently large to prevent the 
tidal breakup of stars passing through the IMBHB, i.e. 
« ~ [(Mi/m)^/^ -I- (M2/m)^/^]r, where r is the radius of 
the star. For r given by Eq. (|3}, one has 



> 0.8 



Ml 
M^ 



1/3 



+ 



M2 

Mo 



1/3' 



m 



11/30 



Ro (15) 



(for the sake of simplicity we assume a circular binary or- 
bit). It follows from Eqs. (fT3|l - (fT5|l that the smaller the 
mass of the star the larger velocity it, in principle, can 
attain; note also that inequality (|15p sets an upper limit 
on the mass of a main-sequence star that can be ejected 
by a shrinking IMBHB. For a main-sequence star of mass 
m = 8 Mq (the minimum mass of single stars produc- 
ing neutron stars) and assuming that AIi — 500 Mq and 
M2 = 300 M0, one has from Eqs. ^ and JTH) that 
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a > 30 R0, {Voo) ~ 1950 kms"^ and 



3 380 kms"^ 



0.625 



\^ 500 Mo 



30 Rg 



-1/2 



(16) 



where = Mi /{Mi + M2). 

The maximum ejection velocity could be somewhat 
higher if one considers encounters involving massive post- 
main-sequence stars, either a blue supergiant or a bare he- 
lium core. In this case, a massive star can approach the 
IMBHB much closer than a main-sequence star. Although 
the blue supergiant star will lose its hydrogen envelope due 
to the tidal stripping, its helium core can pass within sev- 
eral R0 from one of the binary components without being 
disintegrated by the tidal force. To estimate the maximum 
possible velocity attained by the helium core, one should 
consider the second constraint on the binary separation. It 
follows from the requirement that the gravitational radiation 
time-scale of the shrinking IMBHB (Peter 1964), 



icwR 



/ Ml 



Mo 



\ 300 Mo 



\^5OOM0^ 
-1 / X 4 



30 R, 



(17) 



should be larger than the mean collision time, 

-1 



tcoii ~ 2.5 X 10* yr 



107 pc 

Ml +M2 
800 Mq 



rcl 



5 km s~i 

1 / \ -1 

a 



30 R, 



(18) 



that is, a > 15 R©. From Eq. (|16|l . one has that a helium core 
passing through the IMBHB just before the latter merges 
due to the gravitational emission can attain a speed as large 
as - 4 780 km s"^ 



5 N-BODY SIMULATIONS OF 
HYPERVELOCITY STARS 

In this section, we perform numerical simulations of three- 
body scatterings with IMBHs, both single and in binaries, 
in order to obtain the velocity distributions for the ejected 
stars for the mechanisms described in Sects. B^ and B31 The 
simulations are carried out using the sigmaS package, which 
is part of the STARLAlfl software environment (McMillan 

6 Hut 1996; Portegies Zwart et al. 2001b). For a detailed 
description of the set up of the scattering experiments see 
Gualandris et al. (2005). 



5.1 Exchange encounters between binary stars 
and an IMBH 

First, we focus on interactions in which a binary consist- 
ing of two 8M0 main-sequence stars [with radii of 3.5 R©, 
see Eq. Q] encounters a single IMBH of mass in the range 



|http : //www. manybody. org/manybody/ starlab.html | 
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Figure 1. Branching ratio for the outcome of encounters be- 
tween a binary star and a single IMBH as a function of the IMBH 
mass. The two binary components are assumed to be 8 Mq main- 
sequence stars and the semimajor axis is a = 0.1 au. The differ- 
ent outcomes are: merger (crosses), fly-bys (triangles), ionization 
(stars), exchange (squares) and the subset of exchange results 
with a high-velocity escaper {Vac 700kms~^) (bullets). The 
error bars represent the formal (Icr) Poissonian uncertainty of 
the measurement. 



100 — 1 000 M0 . The relative velocity at infinity between 
the IMBH and the centre of mass of the binary is set to 
5kms~'^, in accordance with typical dispersion velocities in 
YMSCs. The semimajor axis is a = 0.1 au to represent the 
case of the tightest binaries. In Fig.[T] we present the proba- 
bility of different outcomes (branching ratios) as a function 
of the IMBH mass Mbh. For each value of the IMBH mass, 
we perform a total of 2 000 scattering experiments, which 
result either in a fly-by, an exchange, or a merger of the 
components of the binary. Ionizations never take place as 
the binary is too hard to be dissociated by the IMBH (Heg- 
gie 1975). Mergers occur in a large fraction (~40 per cent) 
of encounters due to the small orbital period of the binary 
star and mostly involve collisions between the binary com- 
ponents caused by perturbations from the IMBH. Exchange 
interactions occur in about 30 per cent of encounters, with 
a probability increasing with the IMBH mass. Since the bi- 
nary components have equal masses, the probability of ejec- 
tion for the two stars is equal in exchange encounters. Dur- 
ing such encounters, one of the main-sequence stars is cap- 
tured by the IMBH while the other star is ejected, possibly 
with high velocity. These encounters are the relevant ones for 
the production of hypervelocity stars. If the velocity of the 
ejected star exceeds 700 kms~^, we regard the star as a pos- 
sible progenitor of a hyperfast pulsar (indicated with HFP 
in the figure). The figure shows that more massive IMBHs 
are more likely to eject stars with hypervelocities. 

The distributions of velocities at infinity for the escap- 
ing stars are shown in Fig. [2] for three different values of 
the IMBH mass: Mbh = 1OOM0,5OOM0 and lOOOM©. In 
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Figure 2. Velocity distributions at infinity for escaping stars 
in encounters between an equal mass (m = 8 Mq ) binary star 
and a single IMBH for different values of the IMBH mass: 
Mbh = 100 M0 (left-hand panel), 500 Mq (middle panel), 
IOOOM0 (right-hand panel). The binary semimajor axis is a = 
0.1 au. 



order to obtain stars with velocities Voo ^ 700 km s^^, an 
IMBH more massive than a few hundred solar masses is re- 
quired. In the case of a 500 M© IMBH, about 20 per cent of 
all exchanges result in an escape velocity ^ 700kms~^. The 
fraction increases to 40 per cent for the 1 000 M0 IMBH. 
We note that even an IMBH of mass Mbh = 100 Mq can 
occasionally (in about 1 per cent of all exchanges) produce 
an escape velocity ^ 700kms~^. 

In order to derive the probability of obtaining the 
largest possible recoil velocities (see Sect. l4.2| ). we perform 
additional scattering experiments with the following param- 
eters: mi — 40 Mq, 7712 = 8M0 and a — 0.15 au. The ve- 
locity distributions for escapers are shown in Fig.[3]for three 
different values of the IMBH mass: Mbh = 100 Mq , 500 M0 
and 1 000 M0 . The maximum velocity obtained for each set 
of parameters is consistent with the predictions derived in 
Sect. 14.21 The fraction of encounters resulting in velocities 
larger than 700kms~^ increases from about 7 per cent for 
AfsH = 100 M0 to about 70 per cent for Mbh = 1 000 M© . 

The average velocity of escapers scales as a~^^'^, as can 
be seen in Fig.|4]and as is expected from Eq. (|10p . The figure 
shows the average recoil velocity of escapers (solid symbols) 
as a function of the initial binary semimajor axis for three 
diflterent values of the IMBH mass A'/bh = 100 Mq , 500 Mq 
and 1 000 Mq. The empty symbols indicate the velocity Knax 
for which 1 per cent of the encounters have Voo > Knax- The 
average and the maximum velocities increase with the mass 
of the IMBH, as expected from energetic arguments. 

In our systematic study of the effect of the initial semi- 
major axis of the interacting binary, we performed further 
scattering experiments adopting a homogeneous sampling 
in log a. If the distribution of orbital separations in a star 
cluster is flat in log a, like in the case of young star clusters 



Figure 3. Velocity distributions at infinity for escaping stars in 
encounters between a binary consisting of a primary star with 
mass nil = 40 Mq and a secondary star with mass m2 = 8 Mq , 
and a single IMBH of mass Mbh = 100 Mq (left-hand panel), 
500 M0 (middle panel), IOOOM0 (right-hand panel). For this 
case of unequal-mass binaries, we consider as escapers only the 
least massive stars (m2). The binary semimajor axis is a = 
0.15au. 
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Figure 4. Average recoil velocity of escapers as a function of 
the initial binary semimajor axis in the interaction of a binary 
star with IMBHs of different mass: AfBH = 100 Mq (triangles), 
J\/bh = 500 Mq (squares), Mbh = IOOOMq (circles). Solid 
symbols represent the average velocity obtained from a set of 
2000 scattering experiments while the empty symbols indicate 
the velocity Vmax for which 1 per cent of the encounters have 
Voo > Vmax. The errorbars indicate the Icr deviation from the 
mean. For clarity, we only show them for one data set. 
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Figure 5. Velocity distributions at infinity for ejected stars af- 
ter an interaction between a (40, 8) Mq binary and an IMBH of 
mass Mbh = 100 (left-hand panel), 500 (middle panel) 
and IOOOM0 (right-hand panel). These velocity distributions are 
integrated over the entire range of orbital separations for the ini- 
tial binary. Only secondary stars (m2 = 8Mq) are considered as 
escapers in this figure. 



Figure 6. Average recoil velocity of escapers as a function of 
the initial binary semimajor axis in the interaction between a 
binary consisting of a main sequence star with mass mi = 40 Mq 
and a He core with mass m2 = 5 Mq , and a single IMBH of mass 
Mbh = 100 Mq (triangles), 500 Mq (squares), lOOOMg (circles). 
We consider as escapers only the least massive stars, i.e. the He 
cores. Solid symbols represent the average velocity obtained from 
a set of 2 000 scattering experimens while the empty symbols 
indicate the velocity Vmax for which 1 per cent of the encounters 
have Voc > Knax- The errorbars indicate the la deviation from 
the mean. For clarity, we only show them for one data set. 



(Kouwenhoven et al. 2005), we can superpose the results of 
these experiments in order to obtain the total velocity dis- 
tributions of escapers. The resulting distributions for three 
different values of the IMBH mass are presented in Fig.[Sl 
The average and the maximum recoil velocities increase with 
the IMBH mass, as in previous cases. The distributions ap- 
pear much broader than those in Fig.O as a result of the 
sampling in the semimajor axis. 

In order to validate the theoretical predictions in 
Sect. 14.21 we simulated encounters between an IMBH and bi- 
naries consisting of a main-sequence star and a He core. We 
considered main-sequence stars of mass mi = 40 Mq and ra- 
dius ri = 11 Rq and He cores of mass m2 — 5 Mq and radius 
r2 = 0.6 Rq [see Eq. Q] in binaries of semimajor axis in the 
range 0.15 — 1.0 au. Fig.[S]shows the average recoil velocity 
of escapers (solid symbols) as a function of the initial bi- 
nary semimajor axis for three different values of the IMBH 
mass Mbh = 100Mq,500Mq, and IOOOMq. The empty 
symbols indicate the velocity Knax for which 1 per cent of 
the encounters have Voo > Knax- As in the case of encoun- 
ters between binary stars and a single IMBH, the average 
and the maximum velocities increase with the mass of the 
IMBH. Maximum velocities are somewhat higher compared 
to the case of two main-sequence stars (see Fig. 3)). This is 
due to the fact that, for a fixed initial semimajor axis, a He 
star can get much closer to a BH than a main-sequence star. 
The fraction of very close encounters, however, is small and 
the average velocity of escapers is not substantially higher 
than in the case of main-sequence binaries. 



5.2 Encounters between single stars and an 
IMBHB 

Another exciting possibility, discussed in Sect . 14.31 is medi- 
ated by an encounter between a single star and an IMBHB. 
We simulate these encounters by considering a single star 
with a mass of 8Mq and a radius of 3.5 Rq and binary 
components with masses in the range 100 — 500 Mq, both 
in equal- and unequal-mass binaries. The semimajor axis is 
fixed to a = 30 Rq ~ 0.14 au in all cases and the relative 
velocity at infinity between the single star and the centre 
of mass of the IMBHB is set to 5kms~^, as in the previ- 
ous case. We perform 5 000 scattering experiments for each 
set of parameters. Except for a few mergers (about 1-2 per 
cent of all cases) , the vast majority of the encounters result 
in a fiy-by. With the random choice of impact parameter 
adopted in the code (see McMillan & Hut 1996; Portegies 
Zwart 2001a), typical recoil velocities in such large distance 
encounters are modest. The distribution of recoil velocities 
is strongly peaked at small velocities (Voo < 100 kms~^) but 
shows a long tail towards large velocities (Voo > 700 kms^^). 
The ejection velocities of escapers depend sensitively on how 
close the single star approaches the IMBHB during the en- 
counter. High- velocity ejections are realized only if the star 
approaches the binary within a distance comparable to the 
binary separation. The maximum ejection velocity for this 
scenario is achieved when the incoming star passes through 
the binary system, very close to one of the BHs. The dis- 
tance of closest approach is limited only by the tidal radius 
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Table 1. List of BH masses adopted in the scattering experi- 
ments with zero impact parameter, followed by the average ve- 
locity at infinity for escapers, the velocity Vmax for which 1 per 
cent of the encounters have Vac > V'max Eind the percentage of 
encounters for which the escapers achieve a recoil velocity larger 
than 700kms-l. The maximum velocities should be considered 
as lower limits, as the actual value depends on the random sam- 
pling of the initial conditions in the scattering experiments. 



Ml 






Vmax(l per cent) 


# 


(Mo) 


(Mo) 


( km s ^) 


( km s^-"-) 


(per cent) 


100 


100 


820 


2175 


53 


200 


100 


950 


2380 


66 


200 


200 


1130 


3070 


75 


300 


200 


1260 


3230 


80 


300 


300 


1430 


3875 


83 


500 


300 


1610 


4170 


87 


500 


500 


1845 


4800 


91 



in the gravitational field of the black hole, which is roughly 
given by rt = {MBH/m)^^^r. The fraction of encounters re- 
sulting in velocities larger than 700kms~^ increases from 
about 8 per cent for BHs of 100 Mo to 20 per cent for BHs 
of 500 Mo . 

In order to prove the critical dependance of the ejec- 
tion velocity at infinity of escapers on the distance of closest 
approach to the IMBHB, we perform another set of simu- 
lations with zero impact parameter. The values of the BH 
masses adopted in each set of simulations are reported in 
Table[TJ followed by the average velocity at infinity for esca- 
pers, the velocity Knax for which 1 per cent of the encounters 
have Voa > Vmax and the percentage of encounters for which 
the escapers achieve a recoil velocity larger than 700kms~^ 
computed over 5 000 scattering experiments for each set of 
parameters. The results shown in Table[T] show that very 
large velocities can be achieved in close encounters between 
a single star and a massive BH binary. Encounters with im- 
pact parameters close to zero are, none the less, very rare 
and we conclude that interactions between a stellar binary 
and a single massive BH are more likely to result in the 
ejection of an hypervelocity star. 



6 DISCUSSION 

In this paper, we explored the hypothesis that the hyperfast 
pulsars could be the remnants of a symmetric supernova ex- 
plosion of a high-velocity massive star which attained its 
peculiar velocity (similar to that of the pulsar) in the course 
of strong dynamical three- or four-body encounters in the 
cores of YMSCs. We estimated the maximum velocities ob- 
tainable in these encounters and found that they are compa- 
rable to those measured for hypervelocity stars. We therefore 
argue that the origin of hypervelocity stars could be associ- 
ated with dynamical processes in the cores of YMSCs. Such 
star clusters are expected to be much more common than 
hitherto assumed (see Sect[3l). 

Hypervelocity stars could also originate from other star- 
forming galaxies. The position in the sky, the age and the 
observed velocity of the hypervelocity star HE 0437— 5439 
suggests that it was ejected from the Large Magellanic Cloud 



(LMC) rather than from the Milky-way Galaxy (Edelmann 
et al. 2005). Recently, Gualandris & Portegies Zwart (2007) 
investigated this hypothesis by simulating interactions be- 
tween a binary and an IMBH in a young star cluster in the 
LMC and found that the hypervelocity star HE 0437-5439 
could have been accelerated by one of the young star clusters 
in the LMC, but this would require an IMBH of ^ 1 000 M© . 

The cores of globular clusters may also give rise to eject- 
ing stars and millisecond pulsars with high speed, in particu- 
lar, those globular clusters which have a high central density 
and are rich in recycled pulsars. Some globular clusters may 
even harbour an IMBH, which are known to be able to pro- 
duce hypervelocity stars, and those clusters form a natural 
channel for production of solitary millisecond pulsars, whose 
velocities, in principle, can reach very high values. 

In Sect.|4]and Sect.O we show that main-sequence O- 
type stars cannot be ejected from YMSCs with a peculiar ve- 
locity sufficiently high to explain the origin of hyperfast neu- 
tron stars (unless a > 1 000 M© IMBH is present in the clus- 
ter center), but lower mass main-sequence stars or stripped 
helium cores could be accelerated to hypervelocities. We find 
also that dynamical processes in the cores of YMSCs can 
produce stars moving with velocities of ^ 200 — 400 kms~^, 
which are typical of pulsars (e.g. Hobbs et al. 2005) and the 
bound population of halo B-type stars (Ramspeck, Heber & 
Moehler 2001; Brown et al. 2007), and therefore can con- 
tribute to the origin of pulsar velocities in addition to asym- 
metric supernova explosions and disruption of binaries fol- 
lowing supernova explosions). But even in these cases the 
cross-section for producing such high speeds is small and we 
expect that only a fraction of high-velocity runaways have 
been accelerated by this mechanism. Nevertheless, some of 
velocities observed for pulsars may be explained by the pro- 
posed scenario. Even though we present the cross-sections 
for scattering processes discussed in Sect.|3]and Sect.[S] it is 
still hard to estimate the production rate for the proposed 
mechanism. The main uncertainty at the moment is in the 
number and characteristics of the population of young and 
massive star clusters. 

We argue also that some fast moving neutron stars 
could be the descendants of hypervelocity helium cores, 
which due to their short lifetimes ( ^ 10^ yr) can only have 
been ejected quite recently. As a consequence, the neutron 
star should not have had time to travel a long distance from 
its parent star cluster. If we assume that the ~ 2 Myr old 
pulsar PSR B 1508-1-55 is indeed the remnant of a hyper- 
velocity star, then the pulsar separation from the Galactic 
plane of ~ 2.5 kpc implies that its progenitor was a he- 
lium core (cf. Gvaramadze 2007). The hypervelocity main- 
sequence stars of mass > 8 M© ejected at large angles to 
the Galactic plane end their lives at distances > 10 kpc and 
thereby contribute to a population of halo and intergalactic 
supernovae. 

The proper motion vectors of several pulsars located in 
the high-pressure interiors of their associated SNRs show a 
trend (Ng & Romani 2004) toward alignment with the pulsar 
rotation axes (inferred from the symmetry of toroidal nebu- 
lae surrounding the pulsars) . This trend could be understood 
if pulsars receive at birth a kick directed along their rota- 
tion axes. The observed alignment, however, is not perfect. 
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It is clearly pronounced only for the Vela pulsa][f|. A natural 
explanation of the misalignment (if one adopts that all neu- 
tron stars are kicked along their rotation axes) is that the su- 
pernova progenitor star had an arbitrarily oriented peculiar 
velocity, comparable with or larger than the kick velocity re- 
ceived by its descendant at birth. One possibility is that the 
supernova progenitor was a member of a binary system and 
that the pulsar peculiar velocity is in part due to the recoil 
received by the system after the first supernova explosion. 
Another possibility is that the progenitor star was ejected 
from the parent YMSC due to the processes discussed in 
Sect.m Thus, the hyperfast pulsars produced by supernova 
explosion of hypervelocity stars should not show any align- 
ment between their spin axes and proper motion vectors. 
This inference, however, would be difficult to prove since the 
hyperfast pulsars leave rapidly from the confines of their par- 
ent SNRs [on a time-scale of ~ 1.4 x 10'' v^ggg {Esi/nisu)^^'^ 
yr, where uiooo is the pulsar peculiar velocity in units of 
1000 kms"'^, i?5i is the supernova energy in units of 10""' 
erg, and nisM is the number density of the ambient inter- 
stellar medium] . The only known hyperfast neutron star as- 
sociated with a SNR (RX J0822-4300 in PuppisA; Hui & 
Becker 2006; Winkler & Petre 2007) is a radio-quiet object 
(Gaensler, Bock & Stappers 2000) that makes it impossible 
to infer the direction of its rotation axis. One can also in- 
fer the orientation of the pulsar rotation axis using the radio 
polarization measurements (e.g. Deshpande, Ramachandran 
& Radhakrishnan 1999; Johnston et al. 2005). As expected, 
the proper motion of the hyperfast pulsar PSR B1508-f55 

o o 

shows a large misalignment of ~ 71 (or ~ 19 , depending 
on whether or not the linear polarization of the pulsar ra- 
dio emission is parallel or orthogonal to magnetic field; see 
Chatterjee et al. 2005). 
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